All-electrical measurement of crossed Andreev reflection in topological insulators 
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Using a generalized wave matching method we solve the full scattering problem for quantum spin 
Hall insulator (QSHI) - superconductor (SC) - QSHI junctions. We find that for systems narrow 
enough so that the bulk states in the SC part couple both edges, the crossed Andreev reflection 
(CAR) is significant and the electron cotunneling (T) and CAR become spatially separated. We 
study the effectiveness of this separation as a function of the system geometry and the level of doping 
in the SC. Moreover, we show that the spatial separation of both effects allows for an all-electrical 
measurement of CAR and T separately in a 5-terminal setup or by using the spin selection of the 
quantum spin Hall effect in an H-bar structure. 
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Many applications of quantum information require en- 
tanglement of quantum states [TH3] . While entanglement 
was achieved using photons [U [5] , it is an ongoing chal- 
lenge to create entangled electrons in solid state devices. 
S-wave superconductors (SC), which couple spin up and 
spin down electrons to Cooper pairs, could be used to 
provide spin entangled electrons by the inverse crossed 
Andreev reflection (CAR) 0[7], a process which splits a 
Cooper pair in the SC into two spatially separated, but 
spin-entangled, electrons in the normal region by the ap- 
plication of a bias voltage. Hence it is interesting to study 
the properties of the CAR and how its magnitude can be 
controlled. A straightforward way to observe CAR is 
by non-local conductance measurements [51 02 ■ However, 
this method has the drawback that CAR is disguised 
by another non-local process called electron cotunneling 
[TU] , which does not involve Cooper pairs and is therefore 
a parasitic process. More involved experimental setups 
could recently detect Cooper-pair splitting circumventing 
the electron cotunneling processes, including additional 
quantum dots [HJH2] or current noise measurements. 
In this paper we make use of the helicity conservation of 
the edge states in quantum spin Hall insulators (QSHI) 
[TH - HB] to achieve a spatial separation of the CAR from 
all other transport channels. Unlike previous works on 
QSHI-SC-QSHI interfaces [T7HT9] we do not restrict our- 
selves to a phenomenological model of the edge states, 
but solve the full scattering problem within a generalized 
wave matching method [201 OH] ■ This method enables us 
to take into account a finite doping in the superconduct- 
ing region, which is often the byproduct of the proximity 
effect. Exactly this mixing of the electron and hole scat- 
tering channels due to evanescent bulk modes in the SC 
results in the spatially localized CAR signal. Further we 
propose setups which allow for direct, all-electrical mea- 
surements of the CAR process, provided that only one 
spin direction contributes to the transport. The spin fil- 
ter can be realized by contacting individual edges of the 



using the non-equilibrium quantum spin Hall effect |22j . 
c.f. Fig. |l(b) Note, the measurement of CAR as well as 
the spin filter destroy the spin entanglement of the un- 
derlying Cooper pair. However, spin-entangled electrons 
could be produced by inverse CAR processes in a setup 
where both edges are coupled to the same reservoirs. 
Without excluding scattering channels, like it was done 
in p-n junctions [19] . we can achieve CAR of up to about 
50% of the total non-local signal, which can be tuned by 
gating the system. 

We consider junctions between two quantum spin Hall 
insulating leads and a superconductor (SC) of length L, 
as it is shown in Fig. 1(a) together with the respective 
dispersions. The transport in the leads is characterized 
by helical edge states, c.f. the solid lines in the disper- 
sion. When the wire has a finite width W, the linear 
dispersion of the edge states acquires a gap due to the 
overlap of the counter propagating edge states at oppo- 
site edges [53]. This so-called mini-gap is schematically 
shown in the electron dispersion for the left and right lead 
in Fig. 
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In the central region, a superconductor is 
placed on top of the quantum spin Hall insulator, which 
induces a superconducting gap and in general dopes the 
system, which is indicated by the energy shift Ci in the 
central dispersion. As a consequence, the superconduct- 
ing gap opens in the bulk of the system, whose energy 
dispersion is indicated by dashed lines in Fig. 1(a) In 



system separately (5-terminal setup, c.f. Fig. 1(a)) or by 



the leads, we will consider no intrinsic doping C\ 3 = 
for simplicity. We will change the Fermi energy Ep in 
the whole structure which would imitate a back gate in 
the corresponding experimental setup (24) . Since the SC 
acts as a reservoir, its chemical potential is fixed and set 
to ground. 

Let us assume a spin up electron is injected from the left 
lead, like shown in Fig. [T] There are four spin conserv- 
ing scattering mechanisms possible. The electron can be 
reflected (R) or tunnel through the sample (T) as an elec- 
tron. To enter the superconducting condensate it needs 
a partner of opposite spin, which ejects a spin up hole 
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FIG. 1: (a) 5-terminal setup. We consider a 
QSHI-SC-QSHI junction. Electron (hole) edge states 
are indicated by red (blue) lines. Above each region we 
show the corresponding band structure schematically. 
Solid lines indicate edge states, dashed lines are bulk 
solutions. For simplicity only electron states are shown 
in the leads. Holes follow from inversion around e = 0. 
The SC is always grounded, (b) Spin selection via the 
H-bar structure of a QSHI. Driving a current / in the 
upper leg of the structure leads to spin up injection to 
the QSHI-SC-QSHI structure in the lower leg. 



either in the left or right lead. These processes are called 
local (AR) or crossed Andreev reflection (CAR), respec- 
tively. Due to the helicity of the edge states R and CAR 
are only allowed on the opposite edge, while AR and T 
take place on the edge of the incoming electron. Neglect- 
ing the doping due to the superconductor, i.e. C2 = 0, 
Adroguer et al. [IB] found that in such superconducting 
tunneling junctions helicity conservation enforces perfect 
AR, when the width (W) and length (L) become large. 
Here we analyze effects of finite size and doping on the 
SC junction. The former effect introduces coupling of 
counter-propagating edge states and inter-lead tunneling 
due to narrow W and short L respectively. The latter 
effect results in the existence of gapped bulk states in 
the SC region. To tackle this problem, we use a gener- 
alized wave matching approach [20] [21] , which solves the 
full scattering problem within hard wall boundary con- 
ditions (HBC). 

Close to the T point (k — 0), HgTe quantum wells can 
be described by an effective 4x4 Dirac model [15]. Us- 
ing the basis (\E1+) , \H1+) , \E1-) , \H1-)) T , with (El) 
and (HI) designating electron and heavy hole sub-bands 
and ± denoting the Kramers' partners, which we here call 
spin for simplicity, the Hamiltonian reads (k± = k x ±ik y , 



k — \ / k 2 . 
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h(h\) =e(k)I + d a (kn)c 
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with d a (k\\) = (Ak x , — Ak y , Ai(k)) T . a a are the Pauli 
matrices and I is the unit matrix in the subband space. 
Here e(k) = C — Dk 2 is the energy dispersion, M.(k) = 
M — Bk 2 is the k dependent bulk gap and A gives the 
strength of the coupling between electron and heavy hole 
states. C = C2 + Ef is an energy shift due to the dop- 
ing by the superconductor (C2) or due to a back gate 
(Ef). The parameters A, B, D and M depend on the 
quantum well width. We will treat the SC, which couples 
time reversed states, in the usual Bogoliubov-de Genncs 
(BdG) formalism. Since we neglected structure and bulk 
inversion symmetry breaking terms [25], % consists of 
two decoupled blocks so that we can restrict our calcu- 
lations to one spin direction. The influence on structure 
inversion asymmetry breaking terms to the proximity ef- 
fect in helical edge states has been studied in Ref. |26j . 
For the upper block, the BdG equation then reads [27] 
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where we defined Hbac- Eq. ^ is correct in the 
absence of time reversal breaking potentials. ip = 
(uei,+ , uhi.+ , Vex,-, — vhi,-) T , where u and v are elec- 
tron and hole like excitations. A is zero in the leads and 
constant in the central region [55] and e is the excitation 
energy relative to the chemical potential in the SC. 
We will consider transport in ^-direction with HBC in 
y-direction. An analytical solution for HBC exists [23] 
for each region, but the scattering problem couples y- 
modes of different regions so that a simple wave match- 
ing is not possible. Here we use the orthogonal set 
0n(j/) = sm[niry/W] (n S N) to Fourier expand 

the eigenstates in a certain part of the slab 
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ip m (x, y) = exp [ik™x] Xn4>n(y) 
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where m labels the different fc™ values. The x™ are f° ur 
component spinors in the basis used in Eq. ([2]). The fc™ 
and the corresponding %™ can be determined by insert- 
ing the ansatz ^ into the BdG Eq. ^ and rewriting the 
obtained system of equations as an eigenvalue equation 
for the fc™, see [5S] for a detailed discussion. While in 
the SC the electrons and holes are coupled, we can dis- 
tinguish electrons from holes in the leads. 
Having the solutions of the individual regions we can 
build up the scattering states for the QSHI-SC-QSHI 
junctions. In the following we will restrict ourselves to in- 
coming spin up electrons from the left side, i.e. the upper 
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edge state in Fig. 1(a) Since the interfaces break transla- 



tional invariance and since the SC couples electrons and 
holes, the scattering states are combinations of solutions 
Q with all allowed k x , weighted by scattering ampli- 
tudes. The amplitudes are calculated by enforcing con- 
tinuity of the wave function as well as the current at the 
interfaces. The transport coefficients R, T, AR and CAR 
can then be obtained from the absolute values squared of 
the amplitudes connecting in and out-going edge states, 
see e.g. Eqs. (A. 14) and (A. 15) in [29 . These coefficients 



therefore are the probabilities that the incoming electron 
undergoes the corresponding scattering process. All the 
other coefficients couple to evanescent modes and do not 
contribute to transport in the leads. However, the trans- 
port by evanescent modes is crucial in the SC region. 
For the following calculations we have chosen the pa- 
rameters from Ref. [S3]: A = 0.3647 nmeV, B = 
-0.7061 nm 2 eV, D = -0.5315 nm 2 eV, M = 
— 10.09 meV and A = 0.5 meV, which is a realistic value 
for the induced superconductivity by a s-wave SC, e.g. 
for bulk Nb A(Nb) sa 1.45 meV at zero temperature [50] . 
We restrict our analysis to the e = regime. 

In Fig. 2(a) we show the transport coefficients as a 
function of W for a constant L = 100 nm. We find that 
R and CAR approach zero for W > 600 nm, which is 
in agreement with Ref. [H], although C2 = —50 meV 
dopes the system far outside the bulk gap. However, 
we do not find quantized AR, since L is relatively small 
and electrons can tunnel through the junction. Increas- 
ing L restores perfect AR. In contrast to Ref. [H] we 
are interested in finite size effects which dominate at 
W < 600 nm. In particular, we observe large fluctu- 
ations in the signal as well as finite R and CAR. The 
latter reaches its maxima at points, where T is relatively 
small. Therefore it is possible that CAR and T have the 
same strength of about 10% of the total signal in this 
configuration. At these points the charge current in the 
right lead becomes minimal and spin currents dominate 
the transport [3T] . Fig. |2(b)| shows the charge density 
for W — 200 nm. Warm colors show electron like states, 
cold colors symbolize an excess of holes. The two inter- 
faces at —L/2 and L/2 are indicated by vertical black 
lines. Inside the leads we find finite density only along 
the edges. On the left side we have an incoming elec- 
tron, which overlaps with a locally reflected hole, and 
a reflected electron on the other edge. The incoming 
electron is either directly reflected as a hole (AR) or en- 
ters the SC. No direct electron reflection happens at the 
first interface [3T]. Inside the SC the electron is reflected 
several times, which explains the resonant behavior at 
small width. Since the electron solution is evanescent for 
= e < A, the resonant behavior due to multiple reflec- 
tions decays with increasing width. A rough estimate for 
the scale on which the resonant behavior can be observed 
can be obtained by the SC coherence length £0 = ^ v f /A. 
For small energies the quadratic terms play a minor role 
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FIG. 2: (a) The transport coefficients as a function of 
the width for length L = 100 nm, C\ = C3 = and 
C2 — —50 meV. We see that the finite size effects 
increase for smaller width, while at large width CAR 
and R are suppressed, (b) The charge density at width 
200 nm (cf. black line in (a)), i.e. at large CAR. Warm 

(cold) colors and positive (negative) numbers 
correspond to electron (hole) densities. We see that the 
incoming electron edge state (left lead upper edge) is 
scattered through the SC bulk modes into a spatially 
separated cotunneled electron and the CAR hole edge 
state on the right side. 



[3"2"] and the approximation £0 rj A/ A = 729 nm holds 
and all resonant behavior should take place at W < £o- 
In Fig. |2(a)| R and CAR are suppressed earlier, which 
is due to the fact that the particle does not take the di- 
rect way, but traverses the sample several times. In the 
right lead the solutions of T on the upper edge and of 
CAR on the lower edge are spatially separated due to 
helicity conservation. The separation works, as long as 
the edge states, which decay like exp[— Xy] away from the 
edges, do not overlap. This is fulfilled for our setup since 
1/A w 20 nm <C W. This separation allows for a direct 
observation and electrical measurement of the CAR. Be- 
fore we describe the proposed detection schemes, let us 
analyze the occurrence of finite CAR in greater detail. 
In Fig. |3(a) we show the transport coefficients as a func- 
tion of the SC's length L at fixed W = 200 nm. Here 
we observe that for long L the probability of AR goes to 
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FIG. 3: (a) The transport coefficients as a function of L 
for W = 200 nm. In addition to the resonant behavior 
we observe an exponential decay on the length scale of 

£q. The constant c = 0.18 was chosen to match the 
CAR signal, (b) The transport coefficients as a function 
of the global doping Ep at fixed geometry W — 200 nm 
and L = 100 nm. 



one. The other transport coefficients decay due to the 
evanescent nature of the particles. The decay lengths 
can be fitted by the coherence length £ (T), £o/2 (R) 
and Co/3 (CAR). This indicates that the particle has to 
traverse the sample several times in order to contribute 
to R or CAR. The need of additional scattering events 
at the interfaces to couple states at opposite edges shows 
that not the direct overlap of the edge states, which is 
suppressed by the large W, but the existence of evanes- 
cent bulk modes, which span over the whole sample, allow 
for finite R and CAR. At Ci = 0, the absence of weakly 
decaying bulk modes prohibits transitions between the 
different edges which can then only be introduced by the 
direct overlap of edge states. However, this overlap would 
destroy the spatial separation of CAR and T in the right 
lead. 

On top of the evanescent decay we find resonant os- 
cillations of the signal. The latter can be understood 
in the context of Fabry-Perot oscillations. As long as 
L <C £o is fulfilled, the signal oscillates on length scales 
of L = 7r/Re(fc™). However, since many different bulk 



modes contribute, the periodicity can not be quantified 
in general. Nevertheless, instead of changing the geom- 
etry, which is experimentally difficult, the Fabry-Perot 
condition can be fulfilled for a fixed L by varying Re(fc™). 
Since ft™ = k™(Ep) we can use gates to realize maximal 
CAR for further investigation. Since the superconductor 
screens the system, gating can be realized most efficiently 
by back gates. In Fig. |3(b)] we present the behavior when 
the back gate changes the Fermi energy in the whole sys- 
tem. For the energy range shown, the leads are not doped 
outside of the bulk gap, which ensures edge state trans- 
port and spatial separation of T and CAR. The signal 
shows the same kind of oscillations like when we changed 
L. Additionally we observe that the signal breaks down 
around E-p w —7.5 meV. The mini-gap opens around 
this energy in the leads, because the overlap of the edge 
states changes as a function of the energy [25] . Close to 
the mini-gap we find a peak in CAR which is generated 
by the direct overlap of the edge states and hence does 
not allow for spatial separation of the transport signal, 
as mentioned before. 

For the calculations we restricted ourselves to one spin di- 
rection, which is crucial for the spatial separation. Here 
we propose two different ways to detect the separation 
between T and CAR: It is possible to contact the edges 
independently in a 5-terminal setup [33] , like it is shown 
in Fig. 1(a) Grounding all but contact 1, where a volt- 



age Vi, is applied, drives only spin up electrons towards 
the junction. Non-local conductance measurements, i.e. 
dlz/dVx and dl^/dVx, then allow for a direct electrical 
measurement of T and CAR, respectively. In Fig. 1(b) 
we propose to use the non-equilibrium quantum spin Hall 
effect in an H-bar structure for spin selection. Applying 
a current / in the upper leg drives electrons from right 
to left in the upper leg. Due to helicity only spin up 
electrons propagate along the lower edge of the upper 
leg. These spin up electrons are transmitted through the 
bridge to the lower leg, where they are injected to the 
junction and scattered in the R, AR, T or CAR edge 
states. 

In conclusion, we propose that the CAR can be mea- 
sured all-electrically in a QSHI - SC - QSHI junction. We 
presented transport calculations showing that finite size 
effects combined with doping the SC by the proximity 
effect enable a significant CAR, which is spatially sepa- 
rated from all other transport channels. We provided a 
handle to tune the magnitude of the signal by use of a 
back gate for given sample geometries. At points where T 
and CAR have the same magnitudes the non-local charge 
currents cancel and we find a pure spin current in the 
right lead without relying on ferromagnetic elements. We 
also argued, that for actual measurements spin selection 
can be achieved by contacting individual edge states or by 
use of the quantum spin Hall effect. Our results should 
be equally relevant to other 2D topological insulators, 
like graphene [2] and thin films of Bi2Se3 [3~4"] . 
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Details on the method 

Here we follow the detailed discussion of the method 
(see also Ref. |Hj). Putting the Fourier ansatz into the 
BdG equation (see Eq. (|3| and Eq. ^ in the main text re- 
spectively) and multiplying from the left by dycj)^ (y) 
we arrive at 
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where all terms of the Hamiltonian proportional to k y — 
—id y are collected in H ky , all constant terms in % const 
and all with k x (k 2 x ) in H k * (U k *). Additionally we 
introduced x^™ = ^™X™- Defining the vectors x m — 
0C=i, X™= 2 , ■ ■ -) T and X ' m = (X'^v xT=i, ■ • -) T , which 
are built from the 4 component spinors Xn an d X 
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Eq. (A.4) can be written as a matrix equation 
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So doubling the dimension of the system of equations 
by introducing x' m allows us to reduce the problem of 
finding fc™ to a linear eigenvalue equation [55. 36 . The 
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sub-matrices in Eq. (A.5) are 
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Here \r m ^ m \ 2 (\r m ,^ rn \ 2 ) and |i mfi , m | 2 (|i m £,J 2 ) are the 
probabilities that the incoming mode m is reflected into 
the electron in mode (hole in mode m^.) of the left 
lead or transmitted into mode m e R (m 1 ^) of the right lead, 
respectively. Enforcing continuity of the wave functions 
and of the currents [Vfci?BdG(fc)] *&(x, y) allows to cal- 
culate all scattering coefficients. In the setup of Fig. [I] 
in the main text only edge state modes are propagating. 
The evanescent waves in the leads do not enter the scat- 
tering matrix. Let m* D from now on be an edge state 
mode, which propagates in D direction. Then we have: 
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Further we used here D± = (D ± B), K n i,n2 



C ± M. The only term coupling 

i = (0m | k y \4> n2 ). For actual com- Using these transport coefficients, we can express the lo- 
cal and non-local conductance Gij(e) — dli/dVj, i.e. the 
conductance from contact j to i, where i,j — 1,...,4 
corresponds to the contacts in Fig. |l(a) in the main text, 
e = eVi is the voltage applied at contact 1. Assuming 
spatially well separated transport signals, the conduc- 
tance at zero temperature can be approximated by [37j 



different modes is r] niri2 = ( 

putation the Fourier series has to be truncated at a suf- 
ficiently high mode N. 

The solutions fc™ can be characterized as propagating 
(k™ € M) or evanescent (Im(k x ) ^ 0). For real fc™ the 
propagation direction can be determined by the sign of 

rW 

v m = / dyip^x, y) [d kx HvdG{k)] kx ^ km ^ m (x, y). 
Jo 

(A.10) 

Analogously evanescent states with lm(k x ) > 
(lm(k x ) < 0) are decaying to the right (left). While in the 
superconductor the electrons and holes are coupled we 
can distinguish electrons from holes in the leads. In the 
following we will denote right (left) propagating/decaying 
solutions by the index m l R (m l L ), where i = e,h for elec- 
trons and holes in the leads. Equipped with this we can 
write the scattering states in each region. To do so we 
label the k x eigenvalues as: k x [ for the part 1 stretching 
from x = —oo to x = -L/2, k" l 2 for — L/2 < x < L/2 

and k x £ for x > L/2, where i = e,h gives the particle 
character in the leads. Our calculations will be restricted 
to zero temperature, where only incoming electrons (not 
holes) need to be considered. In the case where an elec- 
tron is incoming in mode m the states take the form: 

^i(x,y) =i/j m ,i(x,y) + ^ r m e >m <ip m e tl (x, y) 

ml 

+ ^2r m h tm ip m h A (x,y), (A- 11 ) 

ml 

*2(«) y) =^c mR , m t/m R a{ x , v) + ^2 d mLtm i/j mL>2 (x, y) , 

m R m L 

(A.12) 

t m <* t ,m'>Pm H ,3( x > V) + X! t m R ,m' t Pm R ,3( X > f )' 

(A.13) 



h 

G ai (e) = ^T(c) 



G 31 (e) = --CAR(e) 
h 



(A.16) 
(A.17) 
(A.18) 



Moreover, the knowledge of all scattering amplitudes 
determines the full state ^(x,y) up to normalization. 
^(x,y) is the response to an incoming mode ip m ,i(x,y) 
from the left lead. Therefore it allows us to cal- 
culate the non-equilibrium charge density n{x, y) = 
* t (x,2/)A*(x, y), A = diag[l, 1, -1, -1]. For plotting 
the charge density we use the following normalization: 



f w L L 

J o dy&(x = --,y)*(x = --,y) = l. 



(A.19) 



